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SUMMARY 


A  simple  model  of  D-layer  ionospheric  heating  in  the  presence  of  strong,  high 
frequency  (HF)  radio  waves  is  used  to  predict  the  anomalous,  nonlinear  wave  absorption 
due  to  collisional  and  recombination  effects  induced  by  the  indirect  signal.  It  is  found 
that  little  anomalous  absorption  occurs  until  effective  radiated  power  (ERP)  approaches 
100  dBW;  further  increases  in  power  of  several  dB  beyond  this  100  dBW  threshold  are 
frustrated  by  a  comparable  increase  in  self-induced,  one-way  absorption.  This  trend  of 
increasing  absorption  with  increasing  transmitter  ERP  has  considerable  implications  for 
design  of  communication  or  radar  systems  that  use  ultra-powerful,  high-gain  HF 
transmitters. 


STATEMENT  OF  PROBLEM 


Over  the  last  several  years,  proposals  for  very  powerful  HF  radar  transmitters  have 
involved  radiated  power  levels  that  approach  or  exceed  the  threshold  for  nonlinear 
interactions  to  occur  between  the  radio  wave  and  the  ionospheric  plasma  (ERP  >  90- 
95  dBW). 

In  the  ionospheric  E-  and  F-layers  powerful  HF  waves  produce  a  complex  of  wave- 
medium  interactions  due  to  plasma  diffusion  along  geomagnetic  field  lines  [Bloom,  et 
al.,  1991].  In  the  D-layer  (~60km)  powerful  HF  waves  influence  the  medium  largely 
through  nondiffusive  electron  heating.  This  results  in  (a)  change  in  the  electron-neutral 
collision  frequency  and  (b)  change  in  the  electron  number  density  through  temperature- 
dependence  of  electron  recombination  and  attachment  rates.  The  important  consequence 
for  HF  propagation  is  that  (a)  and  (b)  cause  the  incremental  nondeviative  absorption  of 
the  layer  to  increase.  Thus  the  overall  attenuation  suffered  by  a  plane  wave  (ray) 
traversing  the  D-layer  on  a  two-way  path  may  be  markedly  greater  than  that  predicted  by 
the  linear  theory. 

This  is  the  crucial  question  concerning  HF  radar  transmitters:  does  an  increase  in 
ERP  of,  say,  5  dB  (in  the  vicinity  of  the  nonlinearity  threshold  of  the  D-layer)  lead  to  a 
comparable  increase  in  the  overall  wave  attenuation?  As  we  demonstrate,  the  answer  is  a 
qualified  "yes." 

The  incremental  absorption  is  proportional  to  both  the  electron  collision  frequency 
and  to  the  electron  number  density.  In  the  D-layer,  both  of  the  quantities  are 
proportional  to  electron  temperature  (which  itself  increases  with  the  local  electric  field 
intensity).  Using  simple  models  for  these  dependencies,  ambient  daytime  D-layer 
conditions,  and  an  inverse-distance  field  law,  we  have  computed  the  first-order,  field- 
induced  perturbations  in  electron  temperature,  electron  collision  frequency,  and  electron 
density  along  an  18°  ray  path  at  15  MHz.  (We  have  followed  Ginzburg  [1970]  and 
Gurevich  [1978].) 

The  perturbed  values  are  used  to  recompute  a  "modified"  wave  absorption  along  the 
ray  path.  We  have  done  this  for  four  values  of  transmitter  ERP— 90.  95,  100,  and 
105  dBW.  These  approximate  calculations  indicate  that  an  increase  in  HF  power  of  5  dB 
above  the  100  dBW  level  is  met  by  almost  as  much  in  increased  one-way  absorption 
through  the  wave-medium  interaction  effects.  On  a  two-way  radar  path,  therefore,  the 
advantages  of  the  increased  transmitter  power  are  frustrated  by  the  upward  trend  in  total 
D-layer  absorption.  (See  Fig.  1.) 
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Rgure  1 .  Nonlinear  self-effect  of  powerful  HF  wave  in  the  D-layer. 


UNEAR  THEORY 

In  an  isotropic  medium  (e.g.,  the  lower  ionosphere,  neglecting  the  geomagnetic  field 
at  HF  frequencies,  co »  o)^),  the  complex  refractive  index  is  given  by  the  Appleton 
Formula; 


(A) 


X  iXZ 


(B)  where 


V  =  electron  collision  frequency 

=  Np-IZ{fng  (the  "plasma  frequency") 
N  =  electron  density  (cm‘3) 
e  =  electronic  charge 
nig  =  electron  mass 
(0  =  27^,  /is  the  HF  wave  frequency. 


The  equation  for  a  plane  wave  solution  of  Maxwell’s  equations,  traveling  in  the 
x-direction,  is 
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E  = 


or 


> 


E  =  £oexp[)a]-exp 


0)^ 

j)  tor — Ijix 


where  K  =—-Im{n)  measures  the  decay  of  amplitude  per  unit  distance  and  is  called  the 
c 


incremental  absorption  of  the  medium;  it  is  measured  in  nepers/distance,  where 
1  neper  =  8.7  dB. 

Using  (A)  and  (B)  we  obtain 


ic  = 


to  1  XZ 


1  Nv 


c  24  1  +  2^  2eom^c  4o)2+v2 


So-called  nondeviative  absorption  occurs  in  the  D-layer  where  4=1  and  Nv  is  large. 
In  units  of  dB/km,  we  obtain  the  following  expression  for  incremental,  nondeviative 
absorption  [Knapp  and  Fischer,  1970]; 


(vgf-hicohf 


db/km 


U) 


where  g  and  h  are  correction  factors  that  account  for  the  velocity-dependence  of  electron 
collision  frequency. 

The  total  one-way  absorption  for  a  ray  incident  upon  the  D-layer  with  angle  0  is 

Ja(2)^/cos0  dB 
0 

where  0  is  a  more  complicated  function  of  z  for  a  spherical  earth.  (We  have  used  flat- 
earth  calculations  over  relatively  shon  path  lengths.  See  Fig.  2.) 
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6  =  Incidence  angle 


Transmitter 

Rgure2.  Propagation  geometry. 

The  electric  field  will  be  attenuated  by  the  factor  exp[-/4'<^'^^(z)/8.7]  due  to  non- 
deviative  absorption  along  the  ray  path. 

CALCULATION  OF  NONLINEAR  WAVE  ABSORPTION 

In  order  to  predict  possible  seif-attenuation  of  the  strong  HF  wave,  we  must 
investigate  the  dependence  of  the  quantities  v  and  N  in  Eq.  (1)  on  electron  temperature, 
Te. 

We  take  the  height-dependence  of  the  ambient  electron-neutral  collision  frequency  to 
be  -  1.6  X  with  z  in  km.  An  ambient  daytime  profile  for  the  electron 

density  N{z)  is  shown  in  Fig.  3.  These  ambient  values  are  used  to  calculate  from  Eq.  (1) 
the  cumulative  absorption  A"“'(z)  along  the  given  ray  path.  The  peak  electric  field, 
uncorrected  for  temperature-dependent  wave  absorption,  is  £  =  £oexp[-A*’"'(2)/8.7]  The 
peak  field  amplitude  Eq  is  the  simple,  free-space,  far-field  approximation 


where  PG  is  the  effective  radiated  transmitter  power  multiplied  by  the  antenna  gain;  Zq  is 
the  impedance  of  free  space;  and  R  is  the  distance  along  the  ray  path  from  the  source. 
This  quantity  £(r)  is  subsequently  used  to  compute  a  first-approximation  to  the  change  in 
electron-temperature  at  each  height  z. 

In  the  absence  of  heat-conduction  and  diffusion,  the  electron  temperature  at  a  given 
height  z  is  governed  by  the  differential  equation  [Field  and  Warber,  1985]: 
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Height  (km; 


Rgure  3.  Ambient  daytime  electron  density  profile  Ng{^. 


where 

2  3KTQmg5{aP'  +  ) 

Ep  ■■  ^  , 

Ep  is  the  “characteristic  plasma  field”  [Gurevich,  1978] 

Tg(f)  is  the  electron  temperature 

T  is  the  time  after  the  application  of  the  field  E 

K  is  Bolzmann’s  constant 

7*0  is  the  ambient  electron  temperature 

nig  is  the  electron  massS  is  fractional  energy  lost  by  an  electron  in  one  collision 
1)  is  electron  collision  frequency. 
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Regardless  of  the  form  of  the  collision  frequency  v(T),  for  large  times  t » 
field  E  on  the  order  of  Ep  the  steady  state  will  be  achieved: 

=  .  (2) 
2b  Ej 

which  is  the  factor  by  which  the  field  Eq  increases  the  temperature  Tg  to  a  first 
approximation.  This  is  a  nonlinear  implicit  relation  for  Tg,  insofar  as  £^2  depends  on 

V  =  viTg).  Depending  on  the  predominant  coUisional  processes,  vCTg)  will  be 
proportional  to  some  power  of  Tg. 

In  the  ionospheric  D-layer,  electron  collisions  with  neutral  molecules  dominate  and 

V  s  Vg^(Tg).  Unfortunately,  there  is  in  the  existing  literanire  a  panoply  of  seemingly 

dissimilar  formulas  for  For  example,  Gurevich  [1978]  gives  v^^-T  5/6^  but 

Ginzburg  [1970]  uses  -  T  1/2,  Banks  and  Kockarts  [1973]  give  what  is  now  regarded 
as  the  correct  expression  for 


v,„  =  l.8xlO-‘'N„r, 


1+ 


2 

where  is  the  neutral  density  and  the  correction  term  —!==  accounts  for  collisions  with 
O2  molecules. 

In  the  lower  ionosphere,  where  Tg  ~  250°,  this  expression  is  quite  close  to  that  given 
in  Knapp  and  Fischer  [1970]: 


v,„=2.oxio-'V„r“ 


which  we  adopt  for  its  obvious  simplicity  [Field,  1989]. 

Solving  Eq.  (2)  by  iteration,  we  obtain  the  fractional  change  in  electron  temperature 
T  En 

— ^  =  H — ^  and  recalculate  everywhere  the  temperature-corrected  collision  frequency 


^em  =  ^em 


and  the  temperature-corrected  incremental  absorption 


7 


,  4.6x10^  N{2)v,M 
{%msf  +{03hf 

The  comparison,  then,  of  the  cumulative  absorption, 

z  z 

=Ja(z)£b/ cos d  to  =  j a(2)dz  /  cos  0 

0 

illustrates  the  self-attenuation  of  the  strong  HF  wave  £  in  its  traverse  of  the  D-layer. 

This  comparison  of  the  cumulative,  one-way  absorption  of  a  15  MHz  ray  along  an 
18°  path  is  shown  in  Fig.  4.  Four  power-gain  productS“90,  95,  100  and  105  dBW— are 
compared  for  the  corresponding  effect  on  coUisional  absorption.  Figure  5  shows  the 
corresponding  fractional  increase  /  Tg  in  electron  temperature  as  a  function  of  height. 
The  factor  Tg  /  Tq  is  precisely  the  factor  by  which  the  electron  collision  frequency 
increases. 


Rgure  4.  One-way  total  absorption,  1 5  MHz,  18”  ray,  daytime,  corrected  for  temperature 
dependence  of  electron  collision  frequerK^. 


Rgure  5.  Change  in  electron  temperature.  TJ  T^,  15  MHz.  18“  ray.  daytime. 


INCLUSION  OF  ELECTRON-DENSITY  EFFECTS 

An  additional  increase  in  the  incremental  absorption  will  be  produced  by  changes  in 
the  factor  N  in  Eq.  (1).  The  electron  density  N  will  vary  with  electron  temperature 
through  the  attachment  and  recombination  rates. 

From  Knapp  and  Fischer  [1970]  we  obtain  the  following  expression  for  the  steady 
state  D-layer  electron  density  Ngi 

N,  =  ^^jA^+4qa  -  /^2a 

where 


<7  ~  10^  cm'3  sec*  I  (photoelectron  production  rate) 
a  ~  9  X  10*^7^' I  (recombination  coefficient) 

A~5.4xlo2pV-5/2exp[-54O/r,]+1.3xlO-‘5;,[0]  +  2.5xlO-‘2;j[03] 
p  ~  10  dynes/cm2  (pressure  at  80  km)  . 
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A  is  the  electron-neutral  attachment  rate  n[0]  ~  10^^  ^[03]  ~  10^  (molecular  densities 
taken  from  Banks  and  Kockarts  [19731).  The  values  given  here  are  strictly  empirical 
order-of-magnitude  estimates. 

Figure  6  shows  the  dependence  of  Eq.  (3)  on  electron  temperature  Tg.  The  ambient 
temperature  is  Tq  =  300®  K. 


Te  (°K) 


Rgure  6.  Dependence  of  electron  density  on  temperature,  A/( rj/A/( Tg),  daytime  D-layer, 

Tg  =  300. 

The  result  of  including  the  temperature  driven  change  Ng!  N^  'm  the  electron  density 
in  the  total  absorption  Eq.  (1)  is  shown  in  Fig.  7,  for  the  105  dBW  case  only.  The 
electron  density  increases  at  z  =  60  km  by  about  50  percent  due  to  the  temperature 
dependence  of  the  recombination  and  attachment  rates.  This  affects  the  local  incremental 
absorption  and  hence  the  integrated  cumulative  absorption.  Figure  7,  taken  in 
combination  with  Fig.  4,  indicates  that  a  15  dB  increase  in  source  power  at  the  90  dBW 
level  is  almost  offset  by  the  change  in  two-way  collisional  absorption  caused  by  the 
15  MHz  wave  itself. 
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Absorption  (dB) 


Height  of  ray  point  (km) 

Rgure  7.  Corrected  one-way  absorption,  15  MHz,  18®  ray,  daytime. 
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CONCLUSION 


Figure  7  shows  that,  at  the  105  dBW  ERP  level,  total  one-way  wave  absorption  is 
about  8  dB,  nearly  5  dB  greater  than  that  predicted  by  the  linear  theory;  3  dB  of  this 
increase  can  be  attributed  to  temperature  dependence  of  the  electron  collision  frequency, 
and  the  additional  2dB  is  caused  by  temperature  dependence  of  the  electron-ion 
recombination  processes. 

Noting,  from  Fig.  4,  that  the  total  absorption  at  90  dBW  is  nearly  the  same  as  in  the 
linear  theory,  we  see  that  increasing  transmitter  ERP  by  the  additional  15  dBW  yields 
about  16  dBW  additional  absorption  along  the  two-way  path.  This  is,  we  should 
emphasize,  16  dB  more  absorption  than  what  the  linear  theory  predicts.  Thus,  the  point 
of  diminishing  returns  in  transmitter  ERP  is  reached  at  100  to  105  dBW.  The 
implications  of  this  limitation  need  to  be  further  considered  in  high-power,  HF,  radio  or 
radar  system  design. 
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